Segregation and Phase Inversion in a Simple Granular System. 

Akinori Awazu* 

Department of Pure and Applied Sciences University of Tokyo, 
Komaba 3-8-1, Meguro-ku, Tokyo 153-8902, Japan. 



Abstract 

The segregation and the phase inversion are investigated through a simple granular 
system which consists of only two inelastic hard spheres in a square box with an 
energy source. With the variation of the coefficient of restitution, the mass ratio 
between two spheres or the box size, we show that two types of segregated states and 
crossover between them are realized in such a small simple system. PACS number (s): 

Granular materials exhibit various complex phenomena 1,2 . Examples are segregations 
of particles mixtures with different properties which appear by shaking or stirring them, or 
confined them in a horizontally rotating drum, and so on 3-18 . Recently, many studies for 
several segregated patterns and the crossover between them are reported 10-18 . 

In this paper, instead of carrying out simulations with many particles, we choose a simple 
system consisting of inelastic hard spheres with different masses. Although the system is 
so simple, it will be shown to exhibit segregation of particles. This system is expected as a 
simple model of local processes of highly excited granular mixtures. 

Such small systems with elastic hard discs are recently investigated, and are found to show 
a prototype of solid-liquid phase transition, glass transition, and the transition between the 
static and dynamic friction 19-21 . In the following, we show that this simple system realizes 
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the crossover between two segregated states like the change between different segregated 
patterns obtained in recent studies 10-18 . 

The system under consideration consists of two inelastic hard sphere particles with unit 
radius and different masses which are confined in a two-dimensional square box (Fig. f). 
The left-hand wall is set at the origin of the horizontal axis and is in contact with the energy 
source. All walls are rigid and have a length larger than 4. The interaction between a 
particle and the walls without an energy source occurs only through hard-core collisions. 
We give the position of a light and a heavy particle in the horizontal direction xl and xh 
as the distance from the left wall, and the mass mi = 1 and ran — M > 1. The interaction 
between two particles occur through inelastic hard-core collisions with the the coefficient of 
restitution e. A particle hitting the left-hand wall in contact with the energy source with 
the velocity (vh,v v ) bounces back with the velocity (Vh,V v ) (14 > 0). Here, subscripts h 
and v indicate, respectively, the horizontal and vertical directions. In this paper, we employ 
a heat bath of the temperature T as the energy source because of the simplicity. Then, the 
velocity (Vh, V v ) is chosen randomly from the probability distributions PhiVh) an d P V (V V ) 22 : 
tj i-rr \ rriiVh mV 2 

P h (V h ) = -jr- exp( — ^f) (1) 

^(K) = (|^exp(-^), (2) 

where T is the temperature of the heat bath fixed as unity. (We give the Boltzmann constant 
as 1.) 

In the following, we perform the simulation of the above system with the length of walls 
S is = 4.05, = 4.2, and = 4.5. For this range of S, particles are densely packed and the 
collision between particles occur frequently. Now, in order to characterize the distribution 
of two particles, we define the segregation parameter 5 =< xl — xh >t where <> t means 
the time average. If S ~ 0, the particles are not segregated while the sign of 5 gives the 
average configuration of the two particles. 

Now, we draw the phase diagram of the model against {M, e} (M > 1 and e < 1), 
according to the size of S. Figure 2 shows the diagram for (a) S = 4.05, (b) S = 4.2 and (c) 



S = 4.5. Here, 5 > holds in the region with +, 5 < holds in the region with — , and 5 ~ 
holds in the shadow area drawn by multiple of points. (In this paper, we regard the case 
with — (S — d) x 10~ 3 < 5 < (S — d) x 10~ 3 as 5 ~ where d is the diameter of the particle.) 
Independently of S, each phase diagram has following characteristics. For small M, 5 ~ 
or 5 > holds over all the range of e < 1. When M is increased, the system realizes two 
states with 5 > and 5 < and the crossover between them depends on e. In particular, 
critical value e to realize the crossover between states with 5 < and 5 > is independent 
of M for large limit of M. However, crossover points shift to larger e and larger M with the 
increase of S. This means the change of the packing fraction of particles is also relevant to 
the crossover between 5 < and 5 > states. Moreover, the crossover points for S = 4.05 
form a curve with 'N' shape, while the curve becomes smooth with the increase of S like in 
Fig. 2 (b) and (c). 

Figure 3 (a) shows 5 as a function of e foe M = 1.5, M = 2.0, M = 3.65 and M = 6.0 
with S = 4.05. For large M, 5 takes a minimum value at about e ~ 0.825 and the minimum 
value decreases with the increase of M. In Fig. 3 (b), 5 is plotted as a function of M for 
e = 0.825, e = 0.625 and e = 0.55. Here, S has a peak at M ~ 1.3. With the decrease of e, 
5 for larger M becomes larger than that for M at the peak. 

To expalin the existence of the two segregated states with 8 > and 5 < 0, we consider 
two effects determining the sign of 5, respectively; One effect contributes to the increase of 
5 and the other contributes to the decrease of 5 with the decrease of e. 

First, we study the former effect. The light particle's velocity just after the contact with 
the energy source on the left-hand wall tends to be faster than the heavy particle's because 
these velocities are a decrease function of the mass of each particle. On the other hand, by 
iterations of the collision between two particles, their velocities approach with each other 
because e < 1. This implies that the light particle tends to be located farther from the 
left-hand wall than to the heavy particle. Thus, 5 is expected to increase with the decrease 
of e. The contribution of this effect is expected to increase with 1 — e, and as a rough 
approximation, it is assumed to be proportional to 1 — e. 
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Second, we study the effect by which 8 is decreased with the decrease of e. If e is given 
a value close to 1, the approach of two particles' velocities is slow. Then if M >> 1, the 
light particle moves much faster than the heavy one for most of the time and the collision 
between two particles occurs frequently. Then, the light particle behaves like a potential 
barrier for the motion of the heavy particle. In this case, the heavy particle's motion to 
go accross this potential barrier is important to determine the particle distribution. The 
kinetic energy of the heavy particle with the case of xl < xh is smaller than that of xl > xh 
because, in the former case, the heavy particle can not make a contact with the heat bath 
directly and the energy is supplied only by collisions with the light particle. This means 
that the mean velocity of the heavy particle of the case xl < xh is smaller than that of 
xl > xh, while the ratio between velocities is roughly estimated to be e : 1. Then, the ratio 
between the time required to switch from x L < xh to xl > %h and that to switch from 
xl > xh to xl < xh is given as 1 : e. Here, this effect on 5 is prominent only for the case 
with the large collision frequency which is almost propotional to e. In addition, this collision 
frequency increases with the decrease of S. Then, the contribution from the above effect is 
approximately estimated as 5i ~ (e — l)e/C(S) (C is an increase function of S.). 

By the combination of these effects, 5 for large M is given as 5 ~ A5 + B5i = (l — e)(A — 
^fUe). Here, A and B are given as positive constant values. With adequate A, B, and C 
holding < AC(S)/B < 1, 8 takes a negative value for large e (1 > e > AC(S)/B) and a 
positive value for small e (0 < e < AC(S)/ B). This result also gives that the crossover value 
of e between the state with 8 < and that with 5 > becomes larger with the increase of 
S. If M is small, the contribution of 8i is expected little. Thus, 8 > is realized for small 
M. 

In this paper, the mass segregation and the phase inversion are investigated through a 
system which consists of only two inelastic hard spheres in a square box with heat bath. 
With the variation of the coefficient of restitution, the mass ratio or the box size, two types 
of states with different particle distributions and crossover between them are observed for 
the case with a large mass ratio between two particles. The system we studied here may 
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look too small and simple. However, we expect that this system can describe local dynamics 
of highly excited granular systems, and provide a basis for the understanding of particle 
segregations in granular system consisting of many particles. 

The effects of the gravity, the size differences and the friction should be considered in 
order to investigate the generality of obtained phenomena like the crossover between differ- 
ent segregated states 10-18 . Further analytical study of this system to clarify the presented 
behavior as well as the study of systems with three or more particles are necessary in future. 

The author is grateful to K.Kaneko, M. Mizuguchi, H. Hayakawa, H. Nishimori and M. 
Otsuki for useful discussions. This research was supported in part by Grant-in-Aid for JSPS 
Felows 10376. 
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FIGURES 



Fig. 1. Illustration of the system with two inelastic hard spheres and energy source (left). 

Fig. 2. Phase diagrams of the system for each set (M, e) with (a) S = 4.05, (b) S = 4.2 and (c) 
S = 4.2. + and — mean the sign of S, and 5 ~ in the shadow area. 

Fig. 3. (a) 5 as the function of e for M = 1.5, M = 2.0, M = 3.65 and M = 6.0, and (b) 8 as 
the function of M for e = 0.825, e = 0.625 and e = 0.55 with S = 4.05. 
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